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Abstract 

We study the multi-periodic oscillations in the spontaneous emission rate of an atom in a medium 
with refractive index n sandwiched between two parallel mirrors. The oscillations are not obvious in 
the analytical formula for the rate derived based on Fermi's golden rule but can be extracted using 
Fourier transforms by varying the system scale while holding the configuration. The oscillations 
are interpreted as interferences and correspond to various closed-orbits of the emitted photon going 
away from and returning to the atom. This system provides a rare example that the oscillations 
can be explicitly derived by following the emitted wave until it returns to the emitting atom. We 
demonstrate the summation over a large number of closed-orbits converges to the rate formula of 
golden rule. 
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I. INTRODUCTION 



The spontaneous emission process of atoms in the presence of environments is an active 
research area in recent years. It is well known that environments surrounding the atom 
can greatly modify the lifetime of atomic state such modifications are important in 
applications. From a physical point view, the fundamental problem is to understand the 
general characteristics in the modified spontaneous emission rate. Inspired by the similarity 
between the spontaneous emission rate and atomic absorption spectra, it has been shown 
recently that large scale oscillations in the spontaneous emission rate of atoms near an 



dielectric interface and inside a dielectric slab [2j, |3| are present and that a theory similar to 
closed-orbit theory , l3] can be used to understand such large scale oscillations in the 

spontaneous emission rate. Oscillations in the absorption spectra for atoms and negative 
ions in external electric and magnetic fields are associated with electron closed-orbits going 
away from and returning to the nucleus, the oscillations in the spontaneous emission rate are 
associated with what we call "photon closed-orbits" going away from and returning to the 
emitting atom. Following "closed-orbit theory" , the oscillations in the spontaneous emission 
rate are interpreted as interferences between outgoing emitted electromagnetic wave and 
returning electromagnetic wave traveling along various closed-orbits. 

For an atom near an dielectric interface, there is only one oscillation in the spontaneous 
emission rate corresponding to one Cosed-orbi. g The dynamics of this system is quite 
similar to the photodetachment of a negative ion in the presence of a static electric field [8( . 
For an atom inside a dielectric slab, three oscillations were extracted from the numerically 
calculated spontaneous emission rates, and these three oscillations were shown to correspond 
to three photon closed-orbits Q]. In this article we study the oscillations in the spontaneous 
emission rate of an atom in a medium with refractive index n between two large perfect 
parallel plane mirrors. This system offers several advantages over an atom inside the dielec- 
tric slabs 2]: (1) the emission rate can be derived analytically; (2) because of the complete 
reflections of the mirrors, we are able to extract many more oscillations from the emission 
rate and correlate them with emitted photon closed-orbits; (3) one can follow the emitted 
wave and derive a rate formula similar to the closed-orbit theory. 
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II. THE SPONTANEOUS EMISSION RATE FORMULA AND FOURIER 
TRANSFORMS 



The schematic diagram of the system is shown in Fig. 1. An atom is placed between 
two perfect parallel plane mirrors, and a medium with refractive index n fills up the space 
between the mirrors. Denote the distance between the two mirrors by d. We choose the 
coordinate system such that the origin is in the middle of the two mirrors and the z-axis 
is perpendicular to the mirrors. The emitting atom is placed at z, its distance from the 
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upper (lower) mirror is denoted by d\ and 62- Using Fermi's golden rule and fo 
approach developed for the case when the space between two mirrors is vacuum [9], 
the spontaneous emission rate for this system can be derived. When the transition dipole 
moment dyi is parallel to the mirror planes, the emission rate formula is 

o M -22 
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where W vac is the emission rate in vacuum without mirrors, M is the greatest integer part of 
nk d/ir , k = ojq/c is the wave number of the emitted light in vacuum, z is the coordinate of 
the emitting atom. The derivation of formula (1) is given in Appendix A for completeness. 

To investigate the multi-periodic oscillations in the spontaneous emission rate, we follow 
the approach for the dielectric slabs [2j] and use R = (d 2 — di)/(d 2 + di) to indicate the 
configuration of the system. We also use a variable a to indicate the system size relative to 
a standard one. We take the standard distance between the two mirrors as the wavelength 
of the emitted photon in vacuum, that is, d° = Xq = 2ir/ko, then we have the relationship 
z = aRd°/2 , d\ = a(l — R)d°/2 and d 2 = a(l + R)d°/2. In order to compare with the 
spontaneous emission rate in dielectric slabs [2], we also set the refractive index n = 1.49 
appropriate for the poly methyl methacrylate (PMMA) material between the two mirrors. 
In Fig. 2 we show the spontaneous emission rate as a function a for R = 0, 1/3, and 3/5. 
In Fig. 2(a), R — 0, the emitting atom is in the middle of the mirrors, the emission rate is 
like a saw-tooth; In Fig. 2(b), R = 1/3, the distance between the emitting atom and the 
lower mirror is twice of the distance between the emitting atom and the upper mirror, some 
smaller tooth appeared; In Fig. 2(c), the emitting atom is even closer to the upper mirror, 
the distance between the emitting atom and the lower mirror is four times of the distance 
between the emitting atom and the upper mirror, finer tooth appeared. These figures can 



be compared with the more smooth emission rate in dielectric slabs [2|. 

The oscillations in these figures are not obvious but can be extracted using modified 
Fourier transforms defined by 

W(i) = / [W{a) - W bg ]ae ia ^da, (2) 

J a\ 

where Wb g = nW vac is the non-oscillating background emission rate in the medium without 
mirrors. In Eq.(2) W(oc) represents the dependence of the emission rate as a function of the 
scaling variable a . Assuming 7' is the first peak position of the Fourier transform in the 
range [7 m j n , 7 maa; ], the integration range should satisfy (a 2 ~ «i) > 2-71-/7' and the step size 
Act should satisfy Aa < 0.1 x 27r/7 max to ensure numerical accuracy. In our calculations, 
we took ati = 1, «2 = 16, Aa = 0.01. The calculated Fourier transformations of the rates 
in Fig. 2 as defined in Eq.(2) are presented as solid lines in Fig. 3. The numerical value near 
each peak is the extracted peak position. There are five peaks in Fig. 3(a) when R = 0, 
seven peaks in Fig. 3(b) when R = 1/3, and seven peaks in Fig. 3(c) when R = 3/5. Many 
more peaks are present in Fig. 3 compared to the dielectric slab system Q]. 



III. THE OSCILLATIONS AND CLOSED-ORBITS 

We now provide an explanation of the physical meaning for the peaks and give a quantita- 
tive description for the peak positions and peak heights. Over the last ten years, closed-orbit 
theory has been used successfully to understand the oscillations in the absorption spectra 
of atoms and negative ions in external fields 0, ^, f], Q] • The oscillations in the absorption 
spectra correspond to electron closed-orbits, and they are interpreted as interferences be- 
tween outgoing electron and returning electron. The similarity between the oscillations in 
the spontaneous emission rate and the oscillations in the absorption spectra has inspired the 
recent studies of the oscillations in the spontaneous emission near dielectric interfaces 

aa. 

The oscillations in the spontaneous emission rates correspond to closed-orbits of emitted 
photon and they are interpreted as interference effects between emitted photon going away 
from and returning to the atom following various closed-orbits. This view can be made more 
quantitative for the present system. 

Consider the radiation damping of a dipole antenna [ljj . The dipole moment is de~ luJot , 
where ujq is the frequency of oscillation. The radiation damping rate is Wd = P/U where P 
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is the average radiation power of the antenna, U is the energy of the antenna. The radiation 
damping rate can be written as [3] 

W = ^Im(rf* ■ E) = ^[Im(rf* • E ) + lm(d* ■ E ret )\, (3) 

where we have decomposed the electric field E at the position of dipole antenna into a direct 
part Eq and a returning part E ret . This step for the electric field is similar to a step for the 
electron wave function in the derivation of closed-orbit theory . Using r to denote 

the vector of a point relative to the dipole position, the direct electric field is 

. 1 A 7 1 

E = — {(f x d) x r(-) + [d- 3r(r • d)] X - ^)}e«<*-«-*>. (4) 

The returning electric field E ret depends on the environments or the cavity. If the dipole 
antenna is at distance I from the mirror and is parallel to the mirror, E ret can be calculated 
exactly using an imaging method [133. The returning field at the position of the antenna is 

E , = -—d(— + — —) e «?kl-u t) (5) 

where d is a unit vector in the dipole direction. The damping rate for a dipole antenna in 
a medium with an arbitrary refractive index n near one mirror and parallel to the mirror 
plane can be evaluated as 

3 sm(2nk l) cos(2nk l) sm(2nk l) 
W d = W Q - -W [ — + {2nkQl)2 - {2nhl)2 ), (6) 

I — I 

where Wo = fifrta ■ Eq.(6) is a generalization of the vacuum case n = 1 [13j. When the 
distance between the dipole antenna and the mirror is greater than half of a wavelength, 
the last two terms in the square bracket are small enough and can be neglected. The semi- 
classical approximation for the damping rate of a dipole near one mirror becomes 

W d = Wo - IWo 8 ^, (7) 

where S = 2nkol is the action of the emitted photon going from the atom to the mirror and 
back to the atom. This path is a closed-orbit. For a dipole antenna in a medium with an 
arbitrary refractive index n between two mirrors, the semi-classical formula for the damping 
rate can be generalized to 

W d = Wo - \wo j; MSj + <f>j), (8) 

3 3 



where Sj = nk^Lj is the action of the j-th emitted photon closed-orbit, Lj is the correspond- 
ing geometric length of the orbit, <f>j = —rrijix is a phase correction, rrij counts the number 
of reflections by the two mirrors. We identify W d and W in Eq.(8) with the spontaneous 
rate W and Wt, g for the spontaneous emission rate formula in Eq.(l). Eq.(8) is quite similar 
to the formula in closed-orbit theory. The peak positions in the Fourier transform of Eq. (2) 
are given by the action S® of emitted photon closed-orbits with corresponding peak height 
proportional to g/L -, where L°- is the length of the closed-orbit, g is a degeneracy factor 
counting the number of closed-orbit having the same length. Both S® and L° are evaluated 
for the standard system size corresponding to a = 1. We have the scaling relationship 
Lj(a) = aL°- and Sj(a) = olS®. 

The above predictions for the peak positions and peak heights in the Fourier transforms 
(FT) are marked as crosses in Fig.3. We fixed the constant h using one peak height and 
then calculated other peak height using hg/L®. For all the peaks in Fig.3, the predications 
from Eq. (8) for both the peak positions and the peak heights agree well with the numerically 
extracted values from the rate formula in Eq.(l). For the case R = 0, there are five peaks. 
Each peak corresponds to two closed-orbits. The structures of the closed-orbits are shown 
as inserts next to the peaks in Fig. 3(a). The first peak in Fig. 3(a) corresponds to 
two closed-orbits going away from and returning to the emitting atom after being reflected 
back by the upper or the lower mirror once. These two closed-orbits have the same action 
Sl 2 — 2™ = 9.36; the third (fourth) closed-orbit leaves from and returns to the emitting 
atom after reflected by both mirrors. The action for the third or fourth closed-orbit is 
S® 4 = Ann = 18.72. The fifth (sixth) closed-orbit is reflected three times by the mirrors. 
Their actions are S® 6 = 67m = 28.09. Similar discussions can be made for other peaks in 
Fig.3 (a). The double degeneracy for some of the peaks can be removed when the emitting 
atom is away from the middle point. For example, when R — 1/3, the first peak at 9.35 in 
Fig. 3(a) splits to two peaks at 6.23 and at 12.50 in Fig.3(b). The peak at 6.23 corresponds to 
the closed-orbit produced by the reflection of the upper mirror closer to the emitting atom, 
while the peak at 12.50 is produced by the reflection of the lower mirror that is further away 
from the emitting atom. In Table I, we compare the extracted values and the theoretical 
predictions for the peak positions and peak heights in Fig.3. They agree very well. 

In Fig. 4 we compare the peak positions extracted from Eq.(l) using Fourier transform 
(circles) with the analytic predication (solid lines) as a function of R. 16 peaks can be 
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TABLE I: Peak positions and peak heights extracted from Eq.(l) using Fourier transform (FT) 
are compared with the analytic predications of Eq.(8). 

Peak position Peak height 
peak index FT Eq. (8) FT Eq. (8) 



R=0 



1 s * 


9.35 


9.36 


3.5987 


3. 


.5993 


2 nd 


18.72 


18.72 


1.8175 


1. 


.7977 


3 rd 


28.08 


28.09 


1.1957 


1. 


.1984 




37.45 


37.45 


0.90616 


0. 


.8986 


5 th 


46.81 


46.81 


0.72541 


0. 


.7189 



I s ' 


6.23 


6.24 


2.6839 


2. 


.7009 


2'nd 


12.50 


12.48 


1.3692 


1. 


.3461 


yd 


18.72 


18.72 


1.8102 


1. 


.7977 


4th 


24.95 


24.97 


0.65365 


0. 


.6744 


5 th 


31.22 


31.21 


0.55291 


0. 


.5390 


6 th 


37.45 


37.45 


0.90279 


0. 


.8986 




43.67 


43.69 


0.37436 


0. 


.3853 



R=3/5 



I s ' 


3.73 


3.74 


4.5112 


4. 


.5112 


2'nd 


14.94 


14.98 


1.1261 


1. 


.1263 


3 rd 


18.72 


18.72 


1.7920 


1. 


.7977 




22.50 


22.47 


0.7551 


0. 


.7479 


5 th 


33.66 


33.70 


0.50151 


0. 


.4999 


Qth 


37.44 


37.45 


0.8917 


0. 


.8989 


yth 


41.22 


41.19 


0.40507 


0. 


.4082 



clearly identified in the specified range. The numerical and analytic results for the peak 
positions are in excellent agreement. The large number of peaks identified for this system 

nn 

is in contrast with the previous systems with dielectric interfaces [2|, |3|] where only one or 
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three peaks were clearly identified. The difference in the number of identifiable peaks in 
these systems can be explained by the difference between a dielectric interface and a mirror 
in reflecting emitted photon waves. In deriving Eq.(8), it is clear that the peak height is 
proportional to the wave amplitude when it returns to the position of the dipole. For a 
closed-orbit j, it is being reflected rrij times by mirrors or dielectric interfaces. A mirror 
reflects better than a dielectric interface, consequently the returning wave corresponding to 
the same closed-orbit is stronger in a mirror system and many more peaks can be identified. 

Eq.(8) is a sum over closed-orbits. In Fig.5 we show how the sum over closed-orbits 
converges to the rate formula in Eq.(l). Fig.5(a) demonstrates that 4 closed-orbits already 
gives a rather accurate representation of the rate that requires a large number of terms 
in Eq.(l). When the closed-orbits are increased to 16 as in Fig. 5(b), the oscillations are 
greatly reduced. When we use 100 closed-orbits in Eq.(8) as shown in Fig.5(c), the sum over 
closed-orbits in Eq.(8) accurately represents the golden rule rate in Eq.(l). 



IV. CONCLUSIONS 



In conclusions, we have studied the oscillations in the spontaneous emission of an atom in 
a medium with an arbitrary refractive index n between two large perfect plane mirrors. We 
have compared the golden rule formula in Eq.(l) and the summation formula over closed- 
orbits in Eq.(8). The oscillations in the golden rule formula can be extracted when the 
system scale is varied while the system configuration is fixed. The extracted peak positions 
and peak heights agree well with Eq.(8) based on closed-orbits of emitted photon going 
away and returning to the atom. The oscillations in the spontaneous emission rate are 

Euite similar to the oscillations in the absorption spectra described by closed-orbit theory 
For other systems, in principle, we can follow the propagation of the electric 
field along various closed-orbits until they are back to the emitting atom to calculate the 
radiation damping rate using Eq.(3). This would give us a formula of the rate in terms 
of closed-orbits. Our study suggests the rate formula based on Fermi's golden rule and 
the rate formula involving a sum over closed-orbits provide complementary perspectives on 
spontaneous emission process. 
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APPENDIX A: DERIVATION OF FORMULA IN EQ.(l) 



131 ] . The transition rate is 



The emitting rate is given by Feimi's golden rule 

W = % E I < e ' °l - -P- A ^ W I^K - co eq ) (Al) 

a 

where \e, 0) is the system initial state, in which the atom is excited and there is no photon 
in the field; \q, l a ) is the system final state, in which the atom is in final state and there is 
one photon with frequency u a in the field; q is the charge of the electron of the atom, and m 
is the mass of the electron; p is the momentum of the electron, and A(x') is the quantized 
transverse vector potential, and x! is the position of atom; oj qe is the frequency (E e — E q )/h. 
A(x) can be write as 



A W = E \ ^r—{ A «W a « + ^Hl (A2) 

a. 

where 

V 2 A a (x) + k 2 a A a (x) = 0, k 2 a = ^, (A3) 

V Mx) = 0, (A4) 

n is the refractive index of the dielectric medium, and the mode function are chosen to form 
an orthonormal set: 

d 3 xA* a (x) ■ Ap(x) = 5 a p. (A5) 
Using Eqs. (EE]) and flM]) we find 

W = IT, E wf l , IPe, • Zm**W - <*«,) (A6) 



h 2 ' 2m?enn 2 uj 

a 

where p eq = (e\p\q) . Using qp eq = —imu eq d eq and Eq. (A6) we get 



he n 2 ' to, 

a 

where d eq is the electric dipole matrix element 



2 

= ^ E — 14 • ^W^K - "e q ) (A7) 
a 
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Consider our boundary conditions: the mirrors are in z = —d/2 and z = d/2 , and 
therefore the mode function must vanish at this two planes. According to Ref. [h]] we write 
the mode functions as 

A %i {x) = (^ 2 (k x z)smk 3 (z- (A8) 

A^) = l(^) 1/2 ihz cos k 3 (z -i)- ik 3 k\\ sink 3 {z - ^)]e ifc V (A9) 

where fey is the parallel component of k , and k 3 = j^, (j = 1.2. • • •) . We consider d eq is 
parallel the mirror and set as the angle between k\\ and d eq . Then we get 

^^2- - d k 2 - d 

\A a (x)-d eq \ 2 = — [|(fc|| x z) ■ d eq \ 2 sin 2 k 3 (z --) + -^\k v d eq \ 2 sin 2 k 3 (z - -)} (A10) 

= ^[(sin 2 + §cos 2 0)sin 2 fc 3 (^-^)] (All) 
v K 2 

Using Eq. (|A7|) and flAlip . we get 

TTUJ 2 T 2 r°° r 2w 9\d I 2 h 2 

x sin 2 /^- ^)]— 5(w Q - u; eg ) (A12) 
vr< L 2 2|4| 2 ^ /"~ , fc I 



fte n 2 47r 2 V 

&3 



_ /-co 7,2 

£/ ^,ifc| F (i+g) 



x sin k 3 (z — — 



d n n 2 ko 



S(k ]{ - Jn 2 k 2 - k 2 ) (A13) 



„2 



2 ck ck\\ 

tto; 2 L 2 2|4 g | 2 ^ , 2 d nn 2 k 

tt — ^ — > 7r(H — —i; ^ n ) sin k 3 [z (A14) 

he n 2 A7r 2 V ^ V n 2 k 2 -k 2 + k 2j 3V 2 J ck c V ; 



7T^ 2 g L 2 2|rf eg | 2 n ffio^y-^ | J 2 ^ 2 ^„■^2^' 7r ^ 
2 47r 2 ^ c/c c ^ d 2 n 2 fc; 

3 1 

2W Li «f m 1 d 2 



TJ/ d " Wl , J \ ■ 2/3'^ J l[ \ f\ici\ 

W ™iM 2J 1 + sm (-T- - T) (A16) 



lu 3 \d e I 2 

where W vac = 3 e ^ ^3 is the emission rate in vacuum. 
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FIG. 1: Schematic diagram for emitting atom in a medium with an arbitrary refractive index n 
between two perfect parallel plane mirrors. 
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FIG. 2: Spontaneous emission rate (in unit of rate in vacuum) as a function of a = d/dP , where 
d° is the emitted photon wavelength in vacuum. The three cases in (a), (b), and (c) correspond to 
R = (d>2 — d\)/(d2 + d\) = 0, 1/3, and 3/5, respectively. 
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(c) R=3/5 
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FIG. 3: Fourier transforms of the spontaneous emission rate in Fig. 2 (solid lines) are compared 
with the analytical results (crosses) based on closed-orbits. Closed-orbits corresponding to the 
peaks are shown schematically in the inserts. 
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FIG. 4: Peak positions extracted from Eq.(l) using Fourier transforms (circles) are compared 
with the predictions of closed-orbits in Eq.(8) as a function of R. 16 visible peaks can be clearly 
identified for most of R in this range. 
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FIG. 5: The spontaneous emission rate in Eq.(l) (solid lines) and the representation using a finite 
number of closed-orbits in Eq.(8) (dotted lines), (a) 4 closed-orbits; (b) 16 closed-orbits; (c) 100 
closed-orbits. 
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